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As relative systematic frequency uncertainties in trapped-ion spectroscopy are approaching the low 10“^® 
range, motional frequency shifts account for a considerable fraction of the uncertainty budget. Micromotion, 
a driven motion fundamentally connected to the principle of the Paul trap, is a particular concern in these 
systems. In this article, we experimentally investigate at this level three common methods for minimizing and 
determining the micromotion amplitude. We develop a generalized model for a quantitative application of the 
photon-correlation technique, which is applicable in the commonly encountered regime where the transition 
linewidth is comparable to the rf drive frequency. We show that a fractional frequency uncertainty due to the 
2nd-order Doppler shift below \/S.v/v\ = 1 x 10“^° can be achieved. The quantitative evaluation is verified 
in an interleaved measurement with the conceptually simpler resolved sideband method. If not performed 
deep within the Lamb-Dicke regime, a temperature-dependent offset at the level of 10“^® is observed in 
resolved sideband measurements due to sampling of intrinsic micromotion. By direct comparison with photon- 
correlation measurements, we show that the simple to implement parametric heating method is sensitive to 
micromotion at the level of \/S.v/v\ = 1 x 10“^° as well. 


I. INTRODUCTION 

Laser cooled ions in Paul trapsi are excellent candi¬ 
dates for precision spectroscopy)^^— quantum information 
processing^ and quantum simulation.— High trap 
depths allow trapping times ranging from hours to 
months and the strong confinement permits localization 
to within a few nanometers, allowing for precise control 
of energy shifts due to external fields. With trapped ion 
frequency standards approaching fractional frequency 
uncertainties below 10“^^, motional shifts are currently 
among the most significant contributions— Of partic¬ 
ular concern is the 2nd-order Doppler shift arising from 
micromotion, a periodic motion driven by the confining 
radio-frequency (rf) electric field. Methods for the 
precise determination of the micromotion amplitude are 
therefore required in order to minimize its origins and 
to characterize the residual uncertainty. A similar level 
of control over micromotion is required in experiments 
studying the interaction between cold trapped ions and 
ensembles of neutral atomS)ii“— as its kinetic energy 
might otherwise dominate the dynamics. 

The origins and consequences of micromotion as well 
as methods for its determination have been discussed ex¬ 
tensively by Berkeland et al—^ Various other techniques, 
suitable under different experimental conditions, have 
been investigated since thenJ^^— The decreasing overall 
frequency uncertainty in trapped-ion spectroscopy 
prompts for an increasingly detailed understanding 
of the uncertainty in micromotion detection. In this 
article, we experimentally investigate the reliability of 
three of the most commonly used methods, namely 
sideband spectroscopy, the photon-correlation method, 
and parametric excitation. We show that each is capable 


of achieving a fractional frequency uncertainty due 
to the 2nd-order Doppler shift on the order of 10“^°, 
independent of the mass of the ion. A new model 
is developed for photon-correlation signals, which is 
necessary for a quantitative evaluation under commonly 
encountered experimental conditions. It is verified by 
comparison with resolved sideband measurements. We 
derive an expression for a temperature-dependent offset 
seen in resolved sideband measurements due to intrinsic 
micromotion if the ion is not deep within the Lamb- 
Dicke regime. This contribution becomes relevant at the 
level of 10“^® for our parameters. We also discuss the 
applicability of each method during the operation of an 
optical clock, where the disruption of the interrogation 
needs to be minimized. 

This paper is organized as follows: Section HIl reviews 
ion motion in a linear Paul trap, the origins of micro¬ 
motion and the resulting frequency shifts. Section IIIII 
introduces the three evaluated techniques for minimizing 
excess micromotion and experimentally compares their 
resolution limits and suitability for clock operation. The 
results are summarized in section IlYl 

II. MOTION OF TRAPPED IONS 

In the following, we review the motion of a single ion 
in a Paul trap. For brevity, we restrict the discussion to 
linear traps, as this geometry is used in the experiments 
presented here. In a linear trap, rf confinement is only 
used for the two radial dimensions, while the axial 
trapping potential is provided by a static electric field. 
Our results can be directly extended to Paul traps with 
three-dimensional rf confinement. 
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electric field. Equation o contains motion at two 
frequencies: the secular motion at LOi is an oscillation in 
the harmonic trapping potential with the temperature- 
dependent amplitude ui^i, whereas the motion at is 
referred to as micromotion. It consists of an unavoidable 
“intrinsic” contribution from the periodic displacement 
due to the secular motion (IMM), and the so called 
“excess micromotion” (EMM), which occurs whenever 
there is a residual rf electric field E,f at the trapping 
potential minimum. 


FIG. 1. Left: radial plane of a linear Paul trap with asymmet¬ 
ric rf drive and applied dc voltage Ut to remove the degeneracy 
of radial secular frequencies. Black arrows indicate the direc¬ 
tion of the micromotion, which depends on the rf electric field 
lines at the ion position. Right: the electrode geometry of the 
trap used in the experiments presented here. The compensa¬ 
tion voltages Utc and Uec control the ion position in the radial 
plane, as indicated by the arrows. 


Without loss of generality, we assume that the princi¬ 
pal axes of the radial trapping potential are aligned along 
the directions ui and U 2 , as shown in Eig. [U to obtain 
decoupled equations of motion. They are inhomogeneous 
forms of the Mathieu equation^ 


Ui + (ai + 2qi cos(nrft)) 



QEa.a,i 

m 


( 1 ) 


where Q and m are the charge and mass of the ion, is 
the trap drive frequency, E^^^i is the component of static 
electric field along Ui, and the qi and Ui parameters are 
related to the applied voltages. They are defined as 


qi = -92 = 


2QK^{Ur{ 


and 


( 2 ) 


2uj^ AQn,U, 

m(pni ’ 

with the applied rf and dc voltages and Ut, each 
corrected by geometrical factors k of order 1, the 
distance d of the electrodes from the trap center, and 
the secular frequency of the static axial confinement ujz- 
The purpose of Ut is to break the rotational symmetry in 
a controlled way and thereby define the radial principal 
axes. 

The solution to o, to lowest order in qi and at, is^l 
Ui{t) = (uo^i + ui^i cos{uJit)) (l + ^ cos{n^tt + (4) 


O / 2 

with uJi = + Y ’ (5) 

where oc E^^^i is a displacement of the equilibrium 
position from the rf node due to any residual static 


Under the assumption that the radial confinement 
is determined entirely by the ponderomotive potential 
(tti <C gf), the kinetic energy of the micromotion is equal 
to the radial potential energy#^ This doubles the kinetic 
energy from thermal excitation for each direction with 
rf confinement. Assuming an equal temperature T for 
secular motion in all three dimensions, the total thermal 
kinetic energy in a linear trap is therefore 

= IksT , (6) 

where ks denotes the Boltzmann constant. 


Residual rf electric fields can occur due to a displace¬ 
ment ugy of the potential minimum from the rf node or 
due to a phase shift between the rf electrodesii^ In an 
ideal linear trap, micromotion only occurs in the radial 
plane, as there is no axial component of the rf electric 
field. Real traps, however, break the translational sym¬ 
metry due to segmentation and machining tolerances. 
This can give rise to an axial component of Etf, which 
cannot easily be compensated and should therefore be 
minimized by design. This is of particular concern when 
multiple ions are to be stored with low micromotion.— 


The EMM velocity can be derived by integrating 
the equation of motion 


wfi'emm = QE^t COS(f2,ft) , 

and the corresponding mean kinetic energy is 

(E^r^) = = - (^E^t) 

2 m J 


( 7 ) 


( 8 ) 


The motion of the ion gives rise to a 2nd-order Doppler 
shift of 


Ar 




Q 


2mcVL, 


, 5kBT , ^ 
" 2^ ’ 


where E^^^ denotes the total kinetic energy. In addi¬ 
tion, the ion experiences a 2nd-order Stark shift due to 
the non-vanishing mean-squared electric field (E^), which 
can be derived by averaging the square of the quadrupole 
field 


El{ut(t))^UVtEt,)um = 


qtU?pi 

2Q ' 


it) (10) 
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over the classical ion trajectory (SD for all directions of 
rf confinement. The sign and magnitude of the resulting 
shift 



depend on the differential static polarizability CTs of the 
involved statesJ^ For clock transitions in which the 
Stark effect increases the frequency, the proper choice 
of trap drive frequency allows a cancellation of both 
contributionsj ^^i^^ In species with a low differential static 
polarizability of the clock states, such as In+ and Al+^ 
the Stark shift due to a given rf electric field amplitude 
is typically more than an order of magnitude below the 
corresponding 2nd-order Doppler shift. Although we 
therefore quantify our results in terms of the 2nd-order 
Doppler shift, the Stark shift can be derived from the 
same data, given sufficiently precise knowledge of CTs- 



FIG. 2. Spatial phase dependence of the IMM signal. Top: 
ion undergoing secular motion along the k vector of the prob¬ 
ing laser. The indnced micromotion is assnmed to be parallel 
to the secular motion (thin arrows). Center: while the laser 
phase changes linearly with position, the phase of the micro¬ 
motion jumps by tt as the rf node is crossed. Bottom: if the 
amplitude of the secular motion becomes comparable to 
A/4, contributions from either side of the rf node no longer 
cancel entirely. 


respectively. 

III. DETECTION OF EXCESS MICROMOTION 


In the frame of an ion undergoing periodic motion with 
the velocity v = Focos/Orff -I- ip), the Ist-order Doppler 
shift can be expressed as a laser phase modulation with 
modulation index (3 = k ■ where k is the wave 

vector of the laser. For /3 <C 1, the resulting spectrum 
can be approximated by^ 

E{u}) oc Jo{l3)S{u} - lol) 

+ Ji(/3) {5{u: -lol- - 5{lo -ujl + ^rs)) , 

( 12 ) 


where wl is the frequency of the laser in the laboratory 
reference frame and Ji are the Bessel functions of the 
first kind. The two quantitative methods evaluated in 
this work both determine j3, as it is a direct measure 
of the micromotion amplitude along k. The residual rf 
electric field can be derived from /3 via 




( 13 ) 


A. The resolved sideband method 

The resolved sideband method makes use of a tran¬ 
sition with a linewidth of F <C fltf. In this case, the 
micromotion component parallel to k can be determined 
from a measurement of the relative strengths of carrier 
and sideband transitions in Eq. m- 


Depending on F and the available laser power, the 
Rabi frequencies can be determined either from coherent 
population evolution or from steady-state excitation 
rates. In the latter case, decoherence must be well un¬ 
derstood, e.g. solely due to excited state decay, in order 
to quantitatively determine Sufficiently broad 

transitions allow a direct detection of the steady-state 
fluorescence.— 

According to Eq. (U), in the absence of EMM, the 
ion motion contains only Fourier components at tOi and 
However, the phase modulation of the laser elec¬ 
tric field is described by the factor exp(zA: • u{t)), which 
gives rise to a frequency component at 12^ due to its non¬ 
linearity. The amplitude of this component scales with 
qiul ^ and is therefore temperature-dependent. Figure [2] 
illustrates a more intuitive approach: If it is not deep 
within the Lamb-Dicke regime,— the ion samples a con¬ 
siderable range of laser phases during its secular motion. 
This phase change can prevent the cancellation of the 
out-of-phase IMM contributions from either half of the 
secular period. Geometrical considerations for the two- 
dimensional case are discussed in appendix along with 
a quantum-mechanical derivation which shows that for a 
Fock state with n phonons, the signal can become as big 
as 


Ho 


jvt (2n+ 1) + e>(?7f) , 


(15) 


H±i _ Ji(/3) ^ ^ 
Ho JoiP) 2 


(14) 


where Ho and H±i denote the Rabi frequency of the 
carrier transition and Ist-order micromotion sideband. 


where rji = k^^h/{2muji) is the Lamb-Dicke parameter of 
the mode parallel to k. The signal for a thermal state can 
be derived from this expression as an ensemble average 
of Rabi oscillations. 
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FIG. 3. Illustration of the model used to derive the photon- 
correlation signal. The atomic transition is described as a 
classical damped oscillator (see Eq. 113), the spectral re¬ 
sponse of which is shown by the solid (amplitude) and dashed 
(phase) lines. In its frame of reference, micromotion causes 
an apparent phase modulation of the laser, which produces 
sidebands (vertical bars, see Eq. (1121) ). The quantities r,flrf 
and A are indicated in units of T. 


B. The photon-correlation method 

The photon-correlation method determines the mod¬ 
ulation index /3 using a transition that violates the re¬ 
solved sideband condition T <C It is based on the 
correlation between the trap drive voltage and the fluo¬ 
rescence modulation due to which is observed 

by recording the time delay between each detected pho¬ 
ton and the succeeding rf zero crossing. A histogram of 
the observed delays resembles the time-inverted distribu¬ 
tion of photons within one rf period. It has the form 

S{t) = So + AS cos(r2rft — ip) , (16) 


where Sq is proportional to the mean fluorescence, the 
observation time and the bin size, and the normalized 
modulation amplitude AS/Sq and the phase p are used 
to determine the amplitude and phase of the EMM. 


The method has been described in detail by Berkeland 
et alJ^ for the case when T ^ i.e. assuming that 
the scattering rate is in steady state at every instance of 
the rf cycle. However, with typical rf drive frequencies 
on the order of a few 10 MHz, this assumption cannot be 
made for e.g. the S P transitions in alkali-like ions. 
We derive in the following the signals obtained by the 
photon-correlation method in a way that is also valid in 
the regime T « 


In the limit of low intensity I <C Aat, the atomic transi¬ 
tion can be described as a classical damped harmonic os¬ 
cillator with resonance (angular) frequency wq and damp¬ 
ing rate T. Its frequency response to an excitation at 
frequency uj is^ 


A{oj 


_ 1 - Wo) 

- 2 (co - Wo)2 + (^)2 


(17) 


as shown in Eig. |3l The detected fluorescence is propor¬ 
tional to the power in this oscillator when excited with 
the spectrum (1121) . which in the time domain consists of 
three contributions: 


S{AA) 
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ujo)E{oj)P‘^*doj 


oc J2(/3)|A(A)|2 + j2(/3) (|A(A + + |A(A 

'-V- 

So,“dc component” 



+ 2Jo(/3) Ji(^) |A*(A)A(A + Hrf) - A{A)A*{A - a^)! cos (ad + p) 

' -V-' 

AS, “rf component” 

+ 2J^(^/3')\A(^^ + — f^rf)| cos(2r2j.ft H- 


(18) 


with phase 


(^ = arg(A*(A)A(A + af)-dl(A)A*(A-at)) • (19) 


Here, A = wz, — wq is the detuning of the laser (carrier) 
frequency from the atomic resonance. 

The first two terms in (fTSl) correspond to the offset Sq 
and amplitude AS of the observed signal dnj, whereas 
the term that oscillates at 2r2,f is of order Jf(/3) and 
can be neglected for /3 <C 1. This result reflects that the 
process is analogous to performing frequency-modulation 
spectroscopy on the atomic transition)^ 

Equation m shows that for /3 <C 1, the quantity 


AS/Sq oc Ji(/3)/Jo(/3) is a direct measure of the EMM 
amplitude. The proportionality factor depends on the 
usually well-known parameters T and as well as the 
laser detuning A, as shown in Fig. 01 While our model 
agrees with the one derived by Berkeland et in the 
limit <C r, there is a significant deviation when the 
frequency scales become comparable. This is due to the 
fact that the phase difference in the atomic response to 
carrier and sidebands cannot be neglected, and a first or¬ 
der approximation of the line shape over a range of 2H,f 
fails in this regime (see Fig. IS). It is also apparent that 
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FIG. 4. Photon-correlation signal for a fixed micromotion 
amplitude {j3 = 0.085) and varied laser detuning in the 
limit I -C /sat- For the solid curve, F = 27r x 19.6 MHz and 
rirf = 27r X 25.42MHz (indicated by the vertical line). Ex¬ 
perimental values confirm the expected behavior (see section 
IHIDII . The dashed curve shows the result derived by Berke- 
land et al;/® in the limit Ilrt F. The models agree when 
that assumption is valid. 

the highest sensitivity is no longer attained at a detuning 
of — r/2, but rather close to — flrf. Note that the clear 
maximum at a known detuning provides a simple way to 
experimentally determine the resonance frequency, which 
cannot be measured directly on a cooling transitioni^ Its 
exact amplitude and position depend on the laser inten¬ 
sity if saturation cannot be neglected. This dependence 
is discussed in appendix |51 along with further corrections 
to Eq. (fT51) . 

The phase information obtained with the signal is a 
major advantage of the photon-correlation method: the 
separation of in-phase and out-of-phase components in 
measurements with two nonparallel beams allows a full 
determination of the amplitude and orientation of the 
micromotion within the common plane. The phase infor¬ 
mation also allows a distinction between excess micromo¬ 
tion due to a displacement and due to an rf phase shift, 
since the phase of the former changes by tt when the ion 
is moved across the rf node. 


C. Micromotion minimization using parametric excitation 

The third method we evaluate does not measure mi¬ 
cromotion directly, but detects a displacement uoy from 
the minimum of the ponderomotive potential. When a 




T / K 

FIG. 5. Quantities to consider when choosing the laser cooling 
parameters in parametric excitation measurements: Fluores¬ 
cence (top) and Doppler cooling rate (bottom) as a function 
of temperature for laser detunings of —r/2 (solid lines) and 
—Sr/2 (dashed lines) at an intensity of / = O.l/sat. 


modulation of the rf voltage is applied at 2waec/'ra, with 
n G N, the secular motion is excited parametrically at 
a rate proportional to the displacement Unlike the 
other two methods, parametric excitation is not able 
to measure micromotion due to a phase shift between 
rf electrodes. Besides the simple experimental imple¬ 
mentation, a major advantage of this method is that 
it uses the Doppler shift due to secular motion, rather 
than micromotion, to generate a signal. In principle, 
this allows three-dimensional micromotion compensation 
with a single laser beam, as long as it has projections 
onto all principal axes of the trap. 

The signals observed for this method are changes 
in the amount of fluorescence as the modulation fre¬ 
quency is swept over the resonance. Both increases and 
decreases of the fluorescence have been observed, de¬ 
pending on the choice of parameters.— An increase 
in fluorescence typically occurs for low deviations from 
Doppler temperature, whereas fluorescence decreases 
when the temperature is increased further. In order to 
get a monotonic dependence of the signal with respect 
to EMM, it is therefore necessary to operate in the 
regime of fluorescence increase. The corresponding 
temperature range can be extended by detuning the 
cooling laser further from resonance, as shown in Eig. [5] 
(top), which shows the fluorescence based on the value 
of a Voigt profile at detunings of —r/2 and — 3r/2, 
respectively. Since both the heating and Doppler cooling 
rates depend on temperature, experimental parameters 
must be chosen such that equilibrium is eventually 
reached, to avoid runaway heating. As the sensitivity 
of the method increases with equilibrium temperature, 
there is a trade-off between resolution and robustness, 
which can require several iterations while adjusting 
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parameters as the minimum is approached. 


The Doppler cooling rate can be estimated by integrat¬ 
ing the energy loss due to the momentum imparted by 
absorbed photons^ over one period of secular motion: 


I r 

Ecoal = —hkvo — 

L.t 47r. 


n2ir 


sin(t') 


1-f 


2(A —sin(t')) 
r 


rdt' (20) 


with vq = 


keT 


where k is the projection of k onto the principal axis of 
the mode, and we have assumed I <C /sat- Equation (EOl) 
neglects heating due to the recoil from emitted photons, 
which has negligible influence for temperatures far above 
the Doppler limit. Figure [S] (bottom) shows the nu¬ 
merically calculated values of Eq. dsni) as a function of T. 


The heating rate due to parametric excitation has been 
derived by Savard et al4^ as 

TTLUJ^ 

(Ee.c) = e{t')uUt')dt' , (21) 

where the secular frequency is modulated as u!i(t) = 
+ e(t) and r is short compared to the heating 
rate. To give an example, for the experimental param¬ 
eters described in section IIIID2I and = 100 v/m, the 
initial heating rate at T = 0.5 mK is on the order of 
10 “ 2 i Yg experiments, the fluorescence begins to 

drop for E^t > 100 v/m, indicating that the equilibrium 
temperature exceeds 200 mK at this point. 



FIG. 6. Reduced term scheme of The electric dipole 

transition at 369.5 nm is used for Doppler cooling, state detec¬ 
tion, and photon-correlation measurements. Sideband mea¬ 
surements are performed on the 411 nm electric quadrupole 
transition. For ground state cooling, the ^ 05/2 state is de¬ 
pleted via the 1650 nm transition to the short-lived ^P 3/2 
state. 


that the influence on the fluorescence signal can be ne¬ 
glected. Before each scan, micromotion along x is mea¬ 
sured with photon-correlation measurements using the 
horizontal beams and minimized by adjusting the com¬ 
pensation voltage Utc- Axial micromotion can be inferred 
from the same data and minimized by shifting the ax¬ 
ial ion position. After an initial adjustment, no signif¬ 
icant increase in axial micromotion (within an un¬ 
certainty of 40 V“) been observed over the course of 
four months. 


1. Photon-correlation and sideband methods 


D. Experimental comparison 

All of the methods are tested using a single ion. 

Figure [H] shows the relevant levels and transitions. Light 
at 369.5 nm can be applied from two horizontal directions 
(in the x-z plane as defined in Fig. [U at 25 ° and 155 ° 
with respect to the z axis), and vertically. The collected 
fluorescence is split and imaged onto a photomultiplier 
tube (PMT) and an electron multiplying CCD camera 
(EMCCD). The setup is described in more detail in a 
previous publication,^— We use an rf drive frequency of 
Drf = 27 r X 25.42 MHz. In order to find the resolution lim¬ 
its and to test the agreement of quantitative evaluations, 
interleaved measurements comparing two methods each 
are performed while scanning the compensation voltage 
[Zee (see Fig. HD. This shifts the ion almost entirely (to 
within 2°) along the x direction, inducing micromotion 
along y. In the case of sideband and photon-correlation 
measurements, this component is then detected using a 
vertical beam. An additional horizontal cooling beam is 
applied during vertical photon-correlation measurements 
to prevent axial heating of the ion. Its power is two or¬ 
ders of magnitude below that of the vertical beam, such 


For photon-correlation measurements, the ^Si /2 ^ 
^Pi /2 transition (A = 369.5 nm, P = 27r x 19.6 MHz) 
is used. A time-to-amplitude converter measures the 
delay between each photon detected by the PMT and 
the following zero crossing of the rf voltage. Measured 
intervals are then binned and counted by a multi-channel 
analyzer. The amplitude and phase of the micromotion 
are obtained by fitting Eq. (HID to the resulting signal. 
A conservative integration time of 90s is used for each 
measurement, although much shorter intervals could be 
used without significantly deteriorating the resolution, 
as shown in section IHIEI 

Sideband measurements are performed on the 
^Si /2 ^ 05/2 transition (A = 411 nm, T = 27r x 23 Hz), 

addressed by an ultra-stable laser4^ For state detection, 
the fluorescence on the 369.5 nm transition is observed 
with the EMCCD camera. All measurements are 
performed close to the Doppler limit of 0.5 mK at radial 
secular frequencies around 27r x 470 kHz. A spectroscopy 
pulse of varied duration is applied (200 experimental 
cycles each), and a fit with Rabi oscillations of a 
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FIG. 7. Interleaved measurement using the photon- 
correlation and sideband methods (at 0.5 mK) as the ion is 
shifted radially through the trap center. The error bars of 
the photon correlation data are due to shot noise and the 
systematic uncertainties discussed in appendix The error 
bars of the sideband method assume 10 % fit uncertainty of 
the Rabi frequencies. The function fitted to the sideband data 
takes into account the offset due to IMM discussed in section 

MM 



-1 0 1-10 1 


- w,.d/(27r) / kHz 


thermal ensemble is used to extract the ground state 
Rabi frequencies for both the carrier and micromotion 
sideband. 

Figure [7] shows the result of a comparison 
between photon-correlation and sideband mea¬ 
surements. The fitted slopes of E^f with re¬ 
spect to t/ec are Sp^ = (2360 ± 30) (V/m)/V and 
Sab = (2437 ± 29) (V/m)/V for the photon-correlation 
and sideband method, respectively, showing a residual 
mismatch of less than 4%. The optimum values of 
Uec determined by the fits agree to within 4.2 mV, 
corresponding to a difference in 7?^ of (10 ± 4) v/m. We 
currently have no explanation for this discrepancy, which 
corresponds to a 2nd-order Doppler shift of 3 x 10“^^. 
As explained in appendix [BJ radiation pressure during 
the photon-correlation measurements can be excluded 
as a reason. 

While the photon-correlation result at the minimum 
is compatible with zero micromotion, there is a statis¬ 
tically significant offset of about 37’'*^/m when using the 
sideband method at a temperature of T = 0.5 mK. Mea¬ 
surements with varied Doppler cooling parameters show 
a monotonic dependence of on temperature. When 
the ion is cooled to the motional ground state in the di¬ 
rections of rf confinement, no excitation on the sideband 
is observed. We therefore attribute this signal to IMM, 
as described in section ITlI Al The observed amplitude is 
about one third of the maximum value as predicted by 
Eq. (fTKll . As explained in appendix lAl the signal is ob¬ 
served in our geometry due to the lifted degeneracy of the 


FIG. 8. Parametric excitation signals. Top: fluorescence at 
the parametric resonance as the ion is shifted radially through 
the trap (orange circles; mean and standard deviation of 10 
sweeps each). A parabola is fitted to the six innermost data 
points to determine the position of the minimum. Intermit¬ 
tent photon-correlation measurements are used to quantify 
the amount of excess micromotion (blue diamonds). Bottom: 
examples of signals obtained in different micromotion ampli¬ 
tude regimes. Low amounts of micromotion lead to an in¬ 
crease of fluorescence, whereas strong micromotion produces 
a fluorescence decrease. Labels a to d indicate the corre¬ 
sponding data points in the upper graph. 

radial secular frequencies. The determination of the opti¬ 
mal compensation voltage is not affected, since it can be 
interpolated from measurements at higher displacements, 
where excess micromotion dominates. 


2. Parametric excitation 

The parameters for parametric excitation measure¬ 
ments are chosen as a compromise between sensitivity 
and robustness. As shown in section UlI Cl a higher de¬ 
tuning of the cooling laser from resonance widens the 
temperature range for which both fluorescence and cool¬ 
ing rate increase with increasing temperature. The mea¬ 
surements are performed at a detuning A = —ST/2 and 
intensity ///sat ~ 0.01. The parametric resonance is ex¬ 
cited by an amplitude modulation of the rf voltage, which 
modulates the secular frequency as 

Wsec(/) = WsecO (1 + Wq , 


( 22 ) 


Erf.y\ / (V/m) 





















































with Waec.o = 27r X 500 kHz, and modulation depth 
Too = 0.04. The modulation frequency is swept 
over a range of 27r x 2 kHz around Waec.o during 30 s. 
Figure [5] shows a parametric heating measurement as 
the ion is moved along the x direction using the com¬ 
pensation voltage Hec. Each point is derived from the 
maximum fluorescence of the average of 10 modulation 
frequency sweeps, examples of which are shown in the 
bottom part of the figure. Interleaved photon-correlation 
measurements are performed to quantify the amount 
of micromotion along y. It should be noted that the 
two measurements are sensitive to displacements along 
different axes: while the amount of vertical micromotion 
determined by photon-correlation is proportional to the 
X displacement, the strength of parametric excitation 
depends on the displacement along the respective 
principal axis, which is at about —48 ° from the x axis. 
This dependence has been verified by displacing the ion 
by the same amount in different directions. To ensure 
that the minima in Fig. [5] coincide, the y displacement 
of the ion was minimized beforehand using Ut^ and 
photon-correlation measurements with the horizontal 
beams. The fitted minima agree to within AU^c = 2 mV, 
corresponding to a mismatch of = (4.6 ± 3.7) v/m. 
The sensitivity of the parametric excitation method can 
be increased by lowering the trapping potential in order 
to increase the displacement due to a given residual 
static field 

Parametric excitation could also be used for the de¬ 
tection of axial micromotion. However, the contribution 
of the rf field to the axial trapping potential strongly 
depends on the electrode geometry and is ideally neg¬ 
ligible by design, in particular when compared to the 
static axial potential. Exciting axial motion would there¬ 
fore require a much stronger modulation of the rf voltage. 


3. Resolution limits 

Table |T] summarizes the resolution limit we observe for 
a single measurement with each method and gives an 
overview of other published values. 

In sideband measurements, the resolution is lim¬ 
ited by the lowest observable Rabi frequency on the 
sideband and the available laser power. Currently, 
we can achieve fig = 27r x 140 kHz, limited by the 
output power of the second harmonic generation for the 
411 nm light. The sideband Rabi frequency resolution is 
= 27r X 500 Hz, limited by magnetic field fluctua¬ 
tions on the order of 30 nT between experimental cycles. 
The lowest resolvable modulation index due to these 
technical limitations is /3 = 0.007, which corresponds to 
Erf = 2iv/ni. At a temperature of 0.5 mK however, the 
minimum sideband excitation we observe is limited by 
IMM, as described in section UlI D II The fit uncertainty 


in obtaining from measurements using an ion 

at the Doppler limit is about 10 %. Taking all these 
contributions into account, we experimentally observe 
a minimum of Erf = (50 ± 48) ^/m for an ion at the 
Doppler temperature. In the ground state, the expected 
IMM contribution corresponds to less than I v/„i and 
Rabi frequencies can be determined more precisely. 

The photon correlation method is limited by shot 
noise, as discussed in section IHI El It contributes about 
5V/m to the uncertainty in Erf for an integration time 
of 90s. Taking into account the additional uncertainty 
contributions listed in appendix m we experimentally 
observe a minimum of /3 = (2.7 ± 2.7) x 10“^, which 
corresponds to Erf = (7.1 ± 7.2) v/m. 

For the parametric excitation method, the statistical 
uncertainty of the fluorescence peak is given, as shown in 
Fig. |S1 It is dominated by photon shot noise and corre¬ 
sponds to Erf = 23V/ni. 


E. Applicability for optical clock operation 

During clock operation, varying external electric fields 
need to be compensated at regular intervals. If the 
frequency feedback to the clock laser is interrupted while 
the new compensation voltage values are determined, 
the required uncertainty needs to be achieved as quickly 
as possible42 In the following, we discuss the time 
consumption of micromotion compensation with a single 
ion. In a multi-ion clock, an initial evaluation over the 
full extent of the ion crystal is necessary in order to 
determine the frequency shifts of the individual ions. 
However, for crystals with a small extension compared 
to the distance from the electrodes, it is reasonable to 
approximate fluctuating fields as homogeneous. Drifts in 
EMM can therefore be compensated using measurements 
on a single ion as well. In a segmented trap, the Coulomb 
crystal can meanwhile be stored in a separate segment. 

As shown in section IHI D H the quantitative depen¬ 
dence of micromotion on the radial displacement can be 
mapped precisely, see e.g. Fig. [T] A single measurement 
per dimension is therefore sufficient to determine the 
minimum. A controlled displacement is needed to 
resolve the sign ambiguity when sideband or parametric 
heating measurements are used. For photon-correlation 
measurements, this is not necessary due to the contained 
phase information. 

In a photon-correlation measurement, most of the 
uncertainty contributions are relative and therefore 
negligible due to the fact that the micromotion am¬ 
plitude is already low. In this case the dominating 
contribution, by more than an order of magnitude, is 
the statistical uncertainty of AS". In order to evaluate 
its dependence on integration time, we perform a series 
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method 

species 

A /nm 

n,f /MHz 

UJ /MHz 

/3 x 10^ 

Srf/(V/m) 

Bdc/(V/m) 

1 Ai^D 2/^'| X lO^'^ 

ref. 

sideband 

^«Ca+ 

729 

23.5 

3.4 

4 ± 1 

4.2 ± 1.1 

0.42 ± 0.1 

1.3 ± 0.7 

44 


«Ca+ 

397 

55 

4 

0 ± 60 

0 ± 197 

0 ± 10 

0 ± 530 


sideband (T^Doppler) 

172Yb+ 

411 

25.4 

0.47 

17 ± 16 

50 ± 48 

0.6 ± 0.6 

8.5 ± 16 

this work 

photon-correlation 


194 

8.6 

0.065 




0 ± 20 



»»Sr+ 

422 

23.2 


60 

78 


95 



172Yb+ 

370 

25.4 

0.47 

2.7 ±2.7 

7.1 ± 7.2 

0.09 ± 0.09 

0.17 ±0.35 

this work 

parametric excitation 

«Ca+ 


15 

1.2 


22 

1.22 

90 



'“>Ca+ 


14 

1.5 


0 ± 64 

0 ± 6 

0 ± 860 



172Yb+ 


25.4 

0.47 


0 ± 23 

0 ± 0.3 

0 ± 1.8 

this work 

neutral atom probe 

«^Rb+ 


4.17 

0.35 


0 ± 0.34 

0 ± 0.02 

0 ± 0.057 

22 

high finesse cavity 

i3»Ba+ 

493 

5.3 

1.2 

89 ± 3 

11.1 ± 0.3 

1.73 ± 0.05 

15.1 ± 0.9 


repumper Doppler 

^«Ca+ 

866 

25.8 

3.1 

0 ± 8 

0 ± 12 

0 ± 1 

0 ± 8.6 


ion trajectory vs. C/rf 

172Yb+ 


13.2 

0.133 


0 ± 19 

0 ± 0.09 

0±4.65 



TABLE I. Comparison of the sensitivity achieved in a single measurement with each technique and overview of other published 
micromotion measurements. For the methods not discussed in this article, the lowest published values are shown. Bold figures 
indicate the quantity given in the reference. Note that /3 and |Ai/D 2 /*^| quantify EMM independent of the mass. Factors of 27r 
have been omitted in all 17,5 and cu values for clarity. Note that the sideband value quoted for this work is measured at Doppler 
temperature and therefore limited by the IMM signal described in section [til Al 



integration time / s 


FIG. 9. Uncertainty of the photon-correlation method due 
to shot noise as a function of integration time. The mod¬ 
ulation index /3 is calculated for 17,f = 27r x 25.4 MHz and 
A = 369.5nm. The fitted uncertainty is 1.7 x for a 

count rate of 3600 cps. All other contributions to uncertainty 
are typically more than an order of magnitude lower. 


of photon-correlation measurements with a straylight 
signal. Figure |9] shows the standard deviation of fit 
results for AS of 50 measurements each as a function of 
the integration time t, expressed as the uncertainty in /3 
and ifrt for our experimental parameters. The function 
ap{t) = cr/ 3 (l s)/-\/i hts well to the data, indicating 
that the uncertainty is solely due to photon shot noise. 
Since straylight equally contributes to this uncertainty, 
a high signal to noise ratio is helpful. In our case, 10 s of 
measurement per dimension would therefore be sufficient 
to reduce the respective uncertainty in the 2nd-order 
Doppler shift to 1 x 10“^°. Since the signals can be de¬ 
rived from the fluorescence during Doppler cooling, there 
is no fundamental need to interrupt clock operation at 
all when applying this method. For this purpose, cooling 
beams need to be applied from alternating directions to 
allow three-dimensional compensation. With sufficiently 
low trap heating rates, strongly attenuated auxiliary 
beams can be used to ensure cooling of all motional 
modes, as described for the vertical measurements in 
section IIII PI 


In a sideband measurement without ground state 
cooling, the acquisition of a single data point (200 
cycles) takes at least 2 s, limited by the time required for 
state detection and repumping. Total time consumption 
per dimension is therefore on the order of 10 s. Since 
the optimal compensation voltage is extrapolated from a 
measurement with added EMM, it is affected less by the 
IMM contribution and decoherence than measurements 
with minimized EMM. 

Parametric heating measurements need to be per¬ 
formed within the regime of monotonic dependence of 
the fluorescence peak on the micromotion amplitude if 
the minimum is to be deduced from a single data point. 
The choice of parameters therefore requires a trade-off 
between robustness and resolution. As the fluorescence 
value depends on the equilibrium temperature in the 
presence of heating, laser cooling parameters need to be 
well controlled in order to achieve reproducible values. 
On a timescale of hours, we observe a secular frequency 
instability on the order of 100 Hz, which is comparable 
to the width of the parametric resonance. This requires 
either an active stabilization of [/,f or a scan of the 
modulation frequency. Assuming an averaging time of 
10 s per modulation frequency, total time consumption 
is therefore on the order of 10 to 100s. 


IV. SUMMARY AND CONCLUSIONS 

Precision spectroscopy with trapped ions requires an 
accurate determination of motional frequency shifts. In 
this article, we experimentally compare three commonly 
used methods for the compensation and characterization 
of excess micromotion and show that each of them is 
capable of ensuring a 2nd-order Doppler shift well be- 
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method 

advantages 

limitations 

resolution 

time per point 

sideband 

fast acquisition 

measures ratio fli/no 

—¥ common mode effects cancel 

resolution limited by 

— decoherence 

— laser power 

3 X 10“^“ 
a0 — 0.007 

< 10s 

photon- 

correlation 

contains phase information 

Ej-f direction can be inferred, 
rf phase shift / displacement distinguishable 
applicable continuously during clock operation 

resolution limited by 
— photon shot noise 

3.5 X 10“^^ 
a0 — 0.0027 

10 s to 100 s 
(no dead time) 

parametric 

excitation 

simple implementation 
no 3D laser access required 

no quantitative evaluation 
insensitive to rf phase shift 

2 X 10~^“ 

10 s to 100 s 


TABLE II. Summary of the methods evaluated in this work. Resolutions achieved under our experimental conditions are given 
in terms of the second-order Doppler shift. For the quantitative methods, the corresponding uncertainty in the modulation 
index ap is indicated. 


low 10 independent of the ion mass. Table m sum¬ 
marizes the advantages and limitations of these meth¬ 
ods as well as their experimentally obtained single data 
point resolution and time consumption. A new model 
allows us to quantitatively evaluate photon-correlation 
measurements in the common regime where the tran¬ 
sition linewidth is on the order of the rf drive fre¬ 
quency. The obtained micromotion amplitudes agree 
well with those obtained in resolved sideband measure¬ 
ments. While ground state cooling is not necessary for 
the resolved sideband method, a temperature-dependent 
offset on the order of 10“^® limits the resolution if the 
ion is not deep within the Lamb-Dicke regime. It is 
shown to be due to sampled intrinsic micromotion. The 
optimal compensation voltages determined from inter¬ 
leaved scans agree to within less than 5 mV, which cor¬ 
responds to a residual 2nd-order Doppler shift uncer¬ 
tainty of 3.4 X 10“^^ (photon-correlation/sideband) and 
1 X 10“^^ (photon-correlation/parametric heating). The 
corresponding 2nd-order Stark shift due to excess micro¬ 
motion can be deduced from the same measurements, 
given sufficiently precise knowledge of the differential 
static polarizability. In the case of In+ and A1+, this 
contribution is negligiblei^ 
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Appendix A: Sideband method 

In this appendix, we derive the amount of IMM 
sampled by the sideband method, as given by Eq. (HSJ. 
For simplicity, we first consider only motion in one 
dimension, assuming that the secular motion, the 


resulting micromotion, and the k vector are parallel. 
Furthermore, we neglect the pulsation of the harmonic 
oscillator wavefunctions at 

In the adiabatic approximation of Eq. o, micromotion 
can be described as a position-dependent phase modula¬ 
tion with amplitude kxq/2, where x = xo{a + a^) is the 
position operator, and the rf node defines the origin. A 
laser field detuned by appears in the moving frame of 
reference as 

_ ^^gi(“0-l-firf)tgjfc£gi|fc*9cos(nrft) (Al) 

Expanding the rightmost exponential as a series of side¬ 
bands weighted with the respective Bessel functions 
yields 


E{t) = E^^e^i^o+nr!)t^ikx ^ • 

n— — C)0 ^ ' 

(A2) 

Neglecting all terms detuned by > reduces the ex¬ 
pression to 


E{t) = 


/ kxq\ 


(A3) 


The transition matrix element between internal states |g) 
and |e) in the motional state \n) is 


d ■ E{t) 


9 ,n} = —— { n 


Akx 




/ kxq 


(A4) 

where d denotes the transition dipole element. The over¬ 
lap integral of the motional wavefunction evaluates to 


Akx 




/ kxq\ 


= (2n +l) + 0 {if) , 


(A5) 

with the Lamb-Dicke parameter rj = kxQ. It is straight¬ 
forward to include EMM in the above derivation by 
adding an offset to x in the rightmost exponential in 
m- This results in an imaginary second term in 
(IA5I1 with an amplitude of (n|e*^'^|n),5/2 + 0{q^), as 
expected for the EMM sideband. The 7r/2 relative phase 
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FIG. 10. Illustration of the sampled IMM in resolved side¬ 
band measurements for a typical geometry (not to scale). The 
signal depends on the amount of correlation between secular 
motion (gray lines) and micromotion (black arrows) in the 
radial trajectories of the ion when projected onto k, averaged 
over the duration of the pulse. 


between the terms shows that the two components 
cannot interfere and therefore the overall minimum still 
corresponds to compensated EMM. 

Since the term that produces the IMM signal re¬ 
quires correlated micromotion and secular motion com¬ 
ponents along k, its amplitude is maximized in the ge¬ 
ometry assumed above. In typical setups, optical access 
is such that k must be aligned almost perpendicular to 
the quadrupole field lines, as depicted in Fig. [101 The 
relative amplitude of the signal at Doppler temperature 
can then be calculated from the temporal average of the 
product of secular motion and micromotion in the clas¬ 
sical radial trajectory of the ion, each projected onto k. 
For degenerate principal axes, it vanishes completely due 
to symmetry considerations. If the degeneracy is lifted, 
the amplitude depends on the pulse duration compared 
to the difference in secular frequencies, the amount of 
splitting, and the orientation of the principal axes. 


Appendix B: Photon-correlation method 

This appendix lists corrections to Eq. m and ex¬ 
perimental sources of uncertainty for the excess micro¬ 
motion amplitude. For orientation, example values are 
calculated for the ^S' 1/2 O '^Pi /2 transition in ^'^^Yb+ 
at 369.5 nm, which has a natural linewidth of T = 
2tt X 19.6 MHz. For the trap drive frequency, we assume 
Drf = 27r X 25.4 MHz. 


1. RF pickup 

Since there is usually a high amount of power at 
present in the trap drive electronics, it is difficult to pre¬ 
vent parasitic coupling to nearby signals. As with any 
measurement that relies on the demodulation of a signal 



actual excess micromotion amplitude 


FIG. 11. Influence of rf pickup on the photon-correlation sig¬ 
nal. Depending on the relative phase of the disturbance with 
respect to actual excess micromotion, it can either shift the 
minimum position (dash-dotted line) or prevent the minimum 
amplitude from being measured (dashed line). 

at a given frequency, photon-correlation measurements 
are very sensitive to this kind of disturbance. The influ¬ 
ence on the measurement depends on the relative phase 
between the real signal and the pickup and is shown in 
Fig. HI] for the two limiting cases: If the phase difference 
is exactly 7r/2, the two contributions add up quadrati- 
cally, leading to an overestimation of the amount of resid¬ 
ual micromotion when minimized. Due to the fact that 
the phase of the actual micromotion shifts by tt when the 
rf node is crossed, an in-phase pickup signal will shift the 
observed position of the minimum. If the problem can¬ 
not be solved through shielding of the respective com¬ 
ponents, it can be eliminated in data processing. In our 
setup^l we observe a modulation of the PMT detection 
efficiency with which remains at a constant ampli¬ 
tude and phase for weeks as long as no part of either 
the rf circuit or PMT readout electronics is moved. Its 
typical influence is an error in on the order of 20 v/m. 
This kind of disturbance can be removed from the raw 
data after its relative amplitude and phase have been de¬ 
termined by running the measurement cycle with only 
straylight (if necessary for a longer integration time than 
the actual measurements). 


2. Saturation 

While the treatment in section IHIBI has the advan¬ 
tage of providing an analytic expression, the classical 
approach cannot take saturation into account. Experi¬ 
mentally, it is not desirable to measure at low intensities, 
since longer integration times are necessary to obtain a 
given signal to noise ratio. In order to include saturation 
in our model, we numerically integrate the master equa¬ 
tion of a laser interacting with a two-level atom in the 
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detuning / T 


FIG. 12. Numerically calculated photon-correlation signal 
amplitudes as a function of laser detuning when taking sat¬ 
uration into account. The other relevant parameters are 
flrf = 27r X 25.42 MHz, F = 27r x 19.6 MHz, and 13 = 0.085. 


presence of a sinusoidally varied detuning4^ After a suffi¬ 
cient amount of time to avoid the influence of transients, 
a fit to the temporal evolution of the excited state pop¬ 
ulation extracts the signal amplitude AS/So and phase 
(f. Fig. O shows the frequency dependence of AS/Sq 
for a fixed amount of excess micromotion at various in¬ 
tensities. The low intensity result matches the analytic 
expression m- 

The dependence of signal amplitude on saturation 
gives rise to an experimental uncertainty. As a com¬ 
promise between sensitivity and photon shot noise, we 
use I = 0.4 X Aat for our measurements. A first order 
approximation of the numerical results shows that 


<5 



— X (-0.175) X 6 
So 



(Bl) 


As the laser power is actively stabilized, the biggest con¬ 
tribution to i5(///sat) is the uncertainty in the saturation 
power. At 25%, it gives rise to a 1.8% relative uncer¬ 
tainty in AS/So- 


3. Zeeman shifts 

To avoid the existence of dark states, multiple Zeeman 
components of the transition need to be addressed. If 
their frequencies differ by a considerable fraction of the 
transition linewidth, the influence on the signals shown 
in Fig. [12] cannot be neglected. This can be incorporated 
into the model as an ensemble average. For the measure¬ 
ments in section IIIIDi we apply a bias field of 100 /rT 
along k and use linearly polarized light, such that only 
a± transitions are driven. The line is therefore split into 
two equal components that shift in opposite directions 
by the same amount. Figure [T51 shows the effect on the 



FIG. 13. Influence of line broadening effects on the photon- 
correlation signal. Zeeman broadening assumes two compo¬ 
nents, shifted with opposite signs by the amount indicated 
as Au relative to the unperturbed resonance. For Gaussian 
broadening, Au refers to the Gaussian cr. Top: relative max¬ 
imum value of AS/So compared to the unbroadened case. 
Bottom: detuning for which the sensitivity is maximized. 


photon-correlation signal amplitude and the detuning at 
which the maximum occurs, as a function of this shift. 


4. Other line broadening effects 


The transition is broadened further by the Ist-order 
Doppler shift and laser frequency noise. We approximate 
these effects by an ensemble average over a Gaussian dis¬ 
tribution of detunings. Figure [T3| shows the result as a 
function of the standard deviation of this distribution. 


5. Radiation pressure 


The continuous illumination during the measurement 
produces a net force that shifts the equilibrium position 
of the ion by 


Ax = 


hkr^^ 

mcu^ 


(B2) 


where r^c is the photon scattering rate, and ui quanti¬ 
fies the restoring force of the trapping potential in the 
direction of k. For our experimental parameters and 
a; = 27r X 440 kHz, this amounts to Ax = 2.2 nm. Our 
choice of k ensures that the micromotion induced by this 
displacement is perpendicular to k and therefore does not 
affect compensation. 
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6. Further sources of uncertainty 

As shown in section IIIIBi the signal strength for a 
fixed amount of micromotion depends on laser detuning. 
This is a 2nd-order effect if the measurement is done at 
the sensitivity maximum, and in our case it results in a 
relative uncertainty of 

5 f 4^) ~ 4^ X (-4.4 X 10-3 )MHz-2 X . (B3) 

\ JO J ^ Jo 

After the atomic resonance frequency is determined by 
maximizing AS/Sq for a fixed amount of micromotion, 
we stabilize the laser to a wavemeter^ which in turn 
is calibrated with an ultra-stable reference laser every 
600 s. This limits the absolute frequency uncertainty to 
±2 MHz, which leads to a relative uncertainty contribu¬ 
tion of 1.8%. 

The phase information contained in photon-correlation 
signals allows coordinate transformations of measure¬ 
ments done with multiple beams when the k vectors 
cannot be aligned with the axes of interest. However, 
the phase dependence on laser detuning, shown in the 
lower part of Fig. 01 means that frequency fluctuations 
increase the uncertainty of the outcome. 
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